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Q\ • ABSTRACT 

as : 

r^i ! We show that (1) the newly discovered supernova remnant (SNR), GRO J0852- 

| 4642/RX J0852. 0-4622, was created by a core-collapse supernova of a massive star, 

and (2) the same supernova event which produced the 44 Ti detected by COMPTEL 
■ from this source is probably also responsible for a large fraction of the observed 26 Al 

emission in the Vela region detected by the same instrument. The first conclusion is 

(N ■ 

based on the fact that the remnant is currently expanding too slowly given its young 

age for it to be caused by a Type la supernova. If the current SNR shell expansion 
| speed is greater than 3000 km s , a 15M & Type II supernova with a moderate kinetic 

energy exploding at about 150 pc away is favored. If the SNR expansion speed is lower 
qq ' than 2000 km s _1 , as derived naively from the X-ray data, a much more energetic 

supernova is required to have occurred at ~ 250 pc away in a dense environment at 
^ ■ the edge of the Gum nebula. This progenitor has a preferred ejecta mass of < 10M© 

and therefore, it is probably a Type lb or Type Ic supernova. However, the required 

high ambient density of nu > 100 cm~ 3 in this scenario is difficult to reconcile with 
\ the regional CO data. A combination of our estimates of the age/energetics of the new 

SNR and the almost perfect positional coincidence of the new SNR with the centroid 
. £h ! of the COMPTEL 26 Al emission feature of the Vela region strongly favors a causal 

connection. If confirmed, this will be the first case where both 44 Ti and 26 Al are 

detected from the same young SNR and together they can be used to select preferred 

theoretical core-collapse supernova models. 

Subject headings: gamma rays: observations — supernova remnants — supernovae: 
general 

1. Introduction 

Recently, the COMPTEL instrument aboard the Compton Gamma Ray Observatory 
discovered a new supernova remnant (SNR), GRO J0852-4642, as a strong source of 44 Ti 1.16 
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MeV line emission (Eyudin et al. 1998). The new source is located in the Constellation Vela in the 



direction of I ~ 266° and b ~ —1°. Due to the short lifetime of Ti, r = 90.4 ± 1.3 yr ( jNorman et 



al. 1997| ), the fact that COMPTEL can detect this source at all indicates that it is a young SNR 



at a distance less than 500 pc with an age probably less than 10r ~ 900 yr. 

The COMPTEL 44 Ti source was soon identified with a previously unknown shell-type SNR 
in ROSAT data of the Vela region collected during its 1990 all-sky survey ( [Aschenbach 1998] ). 
Designated as RX J0852. 0-4622, the new SNR is seen within the error circle of the COMPTEL 
source, but only at energies > 1.3 keV. The shell has a radius 8 ~ 1°. Spectral analysis reveals that 
X-rays from the majority of the shell are of thermal origin with a temperature of kT = 2. keV, 
which is significantly greater than that of the surrounding area. This is why the shell can be 
identified only at high energies, and this also indicates a young age. The northern limb of the shell 
appears brighter and has an even higher temperature of 4.Tt '^ keV. This component can also be 
fitted by a power law with photon index a = — 2.6^q |. In addition, the ROSAT data gives an 
upper limit of nu < 0.04((f/500 pc) -1 / 2 cm -3 to the density of the ambient medium in which the 
progenitor of RX J0852.0-4622 exploded ( |Aschenbach 1998Q . 



We show in this Letter that, assuming COMPTEL and ROSAT detected the same source, the 
above observed properties can place rather stringent constraints on the nature of the progenitor 
star of GRO J0852-4642/RX J0852. 0-4622 (hereafter GRO/RX J0852). 

The supernova which produced the observed radioactive 44 Ti would also have synthesized 
about the same amount of radioactive 26 Al which has a much longer lifetime of 1.05 x 10 6 yr 
and whose decay produces a 7-ray photon at 1.8 MeV. Indeed, the COMPTEL instrument has 
detected significant 1.8 MeV emission from the Vela region, whose origin has previously been 
attributed to the Vela SNR ( Dichl et al. 1995| ) but later become less certain ( Diehl et al. 1999 ). 



However, the better coincidence of the centroid of the COMPTEL 1.8 MeV feature with the 
location of GRO/RX J0852 prompts us to investigate if it is GRO/RX J0852 that is responsible 
for most of the observed 26 Al emission. 



2. The Nature of the Progenitor 

The observed 44 Ti line flux is F = (3.8 ± 0.7) x 10~ 5 cm" 2 s" 1 (5.6a). For a nominal 44 Ti 
yield of Y = 5 x 1O -5 M per supernova explosion (SNE, [Wooslcy &: Weaver 1994 ; Woosley fc 



Weaver 1995; Thielemann, Nomoto, Hashimoto. 1996| ), the observed 1.16 MeV flux restricts the 



age of the source to be, 

at a distance of 100 pc. In theoretical models, the yield ranges from a low of 5 x 1O^ 5 M to a 
high of 4 x 1O _3 M0. We calculate the constraint on the age versus distance for the range of the 
theoretical yields, as shown by the solid curves in Fig. [|. 
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Fig. 1. — The constraint on the SNR distance and age by average shock velocity and Ti yield 
from a single SNE. 
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While the thermal nature of the bright northern rim is debatable, the thermal X-ray emission 
from most of the shell is undoubtedly generated from shocks as the supernova blast wave runs into 
the ambient medium. The observed shell temperature of 2.5 keV implies that the shock front is 
currently proceeding at a relatively low velocity of 



The validity of such a direct conversion from the electron temperature to the shock front speed 
requires the electrons to be in equilibrium with ions in the post-shock region. For collisionless 
shocks (as is the case for most young SNRs), this condition may not always hold, if, for example, 
the magnetic coupling between the ions and electrons is not effective enough. In such cases, the ion 
temperature, which reflects the true shock speed, could be greater than the electron temperature 
( Scudder 1995| ). Therefore, v s calculated above is probably a lower limit; the real value could be 



greater by as much factor of 2-3. 

Over the lifetime of the SNR, the mean expansion velocity < v s > may be higher than v s , 
depending on its progenitor type. From arguments we present below, < v s > is probably in the 
range of 2000-5000 km s . For given < v s >, the observed shell radius gives another set of 
constraints on the relationship between age and distance, as shown by the dashed lines in Figure |l|. 
We see that the most probable value for the distance to the new SNR is in the range between 100 
and 300 pc, and for the age is between 600 and 1100 years. At a distance of 250 pc, the SNR has 
a radius of about 4.4 pc. 

Now we discuss what possible progenitors are allowed by the observed and derived properties 
of the SNR. We first assume that the environment in which the progenitor exploded has a uniform 
density nu = 0.05 cm -3 , consistent with the limit derived from the ROSAT data. For an SNE 
in such a medium, there is a universal solution for the early stages of SNR evolution ( [McKee "& 



Truelove 1995| ). Immediately after the explosion, the blast wave expands freely during the so-called 



ejecta-dominated phase. The ejecta maintains a constant expansion velocity with a maximum 
value of 

v e = (^f ) 1/2 = 1-3 x 10 4 El( 2 M~q 2 kms~\ (3) 

where M e ^ & is the mass of the ejecta in units of solar masses and E$i is the kinetic energy of the 
SNE in units of 10 51 ergs. For a Type la supernova (SNIa) of M ej Q ~ 1.4, v e is more than 10,000 
km s _1 . For a Type II supernova (SNII) of M e> Q > 10, v e can be much smaller. Therefore, if the 
progenitor of the new SNR is an SNIa, the expansion must have been significantly slowed down 
which requires an old age. If the progenitor star is more massive than 30Mq, the free-expansion 
velocity of an SNII is close to the observed shock front velocity v s . 

After the blast wave sweeps up about the same amount of mass as the ejecta, the SNR 
evolution enters the Sedov- Taylor phase, during which time the SNR shell velocity decreases with 
time as v oc t~ 2 / 5 . The timescale for the onset of the Sedov- Taylor phase is ( McKee &: Truelove 
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1995), 

t ST ~ 570 E^' 2 M e % 6 ng% yr. (4) 

where n^o.05 = ra_H"/(0.05 cm -3 ). Therefore, an SNII spends much more time in the free-expansion 
phase than an SNIa does. Also, a higher ambient density of nu ~ 1 cm -3 would reduce tsT by a 
factor of 2.7. 

It immediately follows from eq. (||) that if the progenitor of GRO/RX J0852 was an SNIa, 
the SNR would have to be in the Sedov- Taylor phase by now because of the observed small shock 
velocity. However, for a typical SNIa explosion with ~ 1 x 10 51 ergs of kinetic energy, it would 
take about g ^ 

tia ~ t S T (£) 5 2 = 1-5 x 10 5 M-*% n~ H % yr (5) 

for the ejecta to slow down from an initial expansion velocity of 11,000 km s" 1 to the observed 
1300 km s . This age can drop to 10 4 yr if the real v s is a factor of 3 higher than that derived 
from the observed electron temperature by eq. (0), but it is still more than an oder of magnitude 
longer than the data suggest. For an SNIa model to work, the only other alternative is then to 
increase the ambient density by a factor of at least 10,000 to > 500 cm -3 . 

The above argument can be more clearly illustrated by considering the theoretical SNR 
expansion velocity as a function of the SNR age shown in Fig. The constant velocity portion 
of the curves represents the free-expansion (or ejecta-dominated) phase, followed by the v oc t -2 / 5 
Sedov- Taylor phase. We have adopted the kinetic energy for an SNIa to be around 1 x 10 51 ergs 



and for an SNII to be in the range 1 — 2 x 10 51 ergs ( Woosley fc Weaver 1994 , 1995; hereafter 
WW94, WW95). The ejecta mass of an SNIa varies from the standard model of 1.4M Q to special 
He-detonated models of 0.8M & (WW94). The ejecta mass of an SNII is in the range between 
8Mq and 4OM0. The shaded area in Fig. ^ represents the parameter space that is preferred by 
the observations. We have allowed a factor of 3 uncertainty in the current shell expansion speed. 

If we take the default ambient gas density as 0.05 cm -3 , we see from Fig. g that the SNIa 
evolution track is totally inconsistent with the observations. To match the current expansion 
speed, the ambient density has to be higher than 500 cm -3 for the SNIa model to be marginally 
acceptable. In this case, we chose a He-detonation model from WW94 with an ejecta mass of 
I.IMq, which has a 44 Ti yield of ~ 2 x 1O~ 4 M0. The mean expansion velocity for this case is 
about 5100 km s -1 at an age of 800 yr. From Fig. g we see that this model with very high density 
does have a solution with a distance ~ 240 pc. 

For the low ambient density condition, on the other hand, we see from Fig. ^ that (1) there is 
a range of SNII models allowed, and (2) all these models are still in the free-expansion phase, so 
< v s >= v e . The model with the highest allowable free-expansion speed of ~ 3700 km s -1 (upper 
dotted line, Fig. |) has E51 = 1.2 and M e , = 15; it produces 5.7 x 1O -5 M of 44 Ti (WW95). 
Fig. [I] then tells us that the source must be at a distance of about 150 pc and an age of about 700 
yr. Other models with smaller expansion speed have much greater progenitor mass of 25 to 40 
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Fig. 2. — The theoretical SNR shell expansion velocity as a function of the SNR age for SNIa, 
SNIb (including SNIc), and SNII models. The input parameters for the SNIa models are £51 = 1.0, 
M e;0 = 1.4, and nu = 0.05 {dashed line) and £51 = 1.6, M eiQ = 1.1, and nu = 500 {dot-dashed 
line). Those for the SNIb model are £51 = 2.5, M eiQ = 10, and nu = 100 {solid line), and for the 
SNII models {n H = 0.05) are £51 = 1.2, M 6i0 = 15 {upper dotted line), and E 51 = 1.2, M 6i0 = 40 
{lower dotted line). The shaded area represents the region that is allowed by the observations. 
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Mq. The lower dotted line of the SNII range in Fig. || corresponds to E51 = 1.2 and M ei = 40. 
However, this model produces only 2.7 x lCr n M of 44 Ti (WW95); in fact, the 44 Ti yield of all 
the models in WW95 which have v e < 3000 km s _1 is less than 1O _9 M0. 

Therefore, if the current SNR expansion speed is less than 2500 km s , the remnant of a 
viable core-collapse supernova model has also to be in the Sedov- Taylor phase which requires a 
large ambient density of nn > 100 cm -3 . For reasons we discuss below, this condition dictates the 
source distance to be at least 250 pc. From Fig. |l| we see that at such a distance the SNR mean 
expansion speed has to be > 5000 km s _1 because the 44 Ti yield of these models is < 3 x 1O~ 4 M 
(WW95). Therefore, we need small ejecta mass and large kinetic energy. The lowest ejecta 
(9.7M0) model of WW95 has a 44 Ti yield of 6 x 10 -5 M© and also a relatively low kinetic energy, 
1.3 x 10 51 ergs. It will produce a mean expansion speed of only 3700 km s _1 for n# = 300 cm -3 . 
A more energetic explosion usually requires a much more massive progenitor and thus even lower 
mean expansion speed. It is possible, however, that a massive star may undergo a supernova 
explosion (resulting in a Type lb or Type Ic supernova, for example) after it loses a significant 
fraction of its envelop because of strong stellar winds. While WW95's calculations did not take 
into account of this effect, we assume that it is conceivable that some core-collapse supernovae 
could have a kinetic energy as much as 2.5 x 10 51 ergs and at the same time an ejecta mass 
of WMq. If so, Fig. [2] shows that such a model (marked as SNIb) can produce a marginally 
acceptable evolution track in an ambient density of 100 cm -3 . The mean expansion speed in this 
case is 5100 km s _1 , consistent with the required distance of 250 pc (Fig. ||). 



3. Discussions 

We have seen that whether the progenitor of GRO/RX J0852 is an SNIa or an SNII, if the 
current SNR expansion speed, v s , is less than 2500 km s , it is almost exclusively required that 
the SNE occurred in a high density region, njj > 100 cm~ 3 for an SNII and > 500 cm~ 3 for an 
SNIa. The distance to the source in either case is about 250 pc and most of the space along the 
line of sight still has a very low density of << 1 cm -3 . We now discuss whether the required high 
density medium exists at the required direction and distance. 

Along the general direction of GRO/RX J0852, the only possible candidate of a high density 
region within 500 pc is the Gum Nebula. While most of the Gum Nebula is filled with ionized, low 
density gas ( [Sahu &: Sahu 1993 ; Reynolds 1976 ), recent reanalysis of the historical HI 21 cm data 



revealed two distinct structures associated with the Gum Nebula ( Reynoso fc Dubner 19971) . One 
is a large, 1.4 x 1O 5 M , disk of neutral gas which is 150 pc in radius and ~ 50 pc in thickness. 
The other is a smaller, thick HI shell of 2.1 x 10 4 M Q with a radius of ~ 60 pc ( pubner et al. 



1992 ), roughly at the center of the disk. Both the disk and the shell share a distance of 500 ± 100 
pc. Therefore, the minimum distance to the edge of the disk is about 250 pc. The mean number 
density of the disk and the shell is 1.6 cm -3 and 0.95 cm -3 , respectively, far below the required 
value for the model progenitor of GRO/RX J0852. Thus, both the SNIa and SNIb models require 
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the presence of additional, much denser gas components at the edge of the HI disk. Because the 
Gum nebula was probably created by some ancient supernova events, it is conceivable that the 
highest density regions are located near the edge of the structure. 



A recent CO survey has shown complex molecular cloud structure in the Vela region (Murphy 



fc May 1991]). The map shows that GRO/RX J0852 falls on the edge of one of the CO hot spots 



but there is no morphological evidence that the cloud is disturbed by the SNR. The distances to 
the components of this so-called Vela molecular ridge are estimated to be 1-2 kpc based on the 
radial velocity profile, further suggesting that they may be unrelated to the SNR. If for any reason 
that these distance estimates were inaccurate and the CO hot spot near GRO/RX J0852 had a 
distance of 250 pc instead of 1-2 kpc, the maximum mean number density of this component would 
be even smaller than the observed 16 cm' 3 ( Murphy Sz May 1991| ), and fall short of the required 



minimum density of 100 cm -3 . Therefore, an extreme density of 500 cm" 3 , along with it the SNIa 
origin of the new SNR, can be ruled out by the HI and CO data; even the SNIb progenitor also 
becomes a difficult proposition. 

If GRO/RX J0852 was created by a SNII event about 800 yr ago, an unavoidable consequence 
is that the same explosion which produced the observed 44 Ti will also synthesize certain amount 
of 26 Al. Specifically, for the SNIb event favored by a small v s the 26 Al yield is 3.6 x 1O~ 4 M with 
a 1.8 MeV line flux of ~ 6.7 x 10 -5 photons s _1 cm -2 at a distance of 250 pc. This flux is about a 



factor of 3 greater than the observed flux observed from the Vela region by COMPTEL flDiehl et 



al. 1995 ). The 15M SNII model favored by the condition v s > 3500 km s , on the other hand, 
has an 26 Al yield of 4.3 x 1O~ 5 M , which would produce a 1.8 MeV line flux of 2.2 x 10 -5 photons 
s _1 cm" 2 , comparable to the observed value. It is thus intriguing to postulate that GRO/RX 
J0852 is also a major contributor to the Vela 26 Al feature. Indeed, recent reexamination of the 
Vela data has argued against previous attribution of the Vela 26 Al flux to the Vela SNR (e.g., 



Diehl et al. 1999|). The centroid of the Vela 26 Al feature, at / = 267° and b = -1°, is ~ 4° away 



from the center of the Vela SNR but falls right on top of GRO/RX J0852. The combination of 
the positional coincidence and the good agreement between the expected 1.8 MeV flux and the 
observed value strongly favors this interpretation. Since the 1.8 MeV Vela feature is probably 
slightly extended, it seems that both the Vela SNR, or a local concentration of the massive-star 



formation tracers ( Diehl et al. 1999 ), and GRO/RX J0852 have contributed while the majority of 



the flux comes from the latter. 



4. Conclusions 

We have investigated the parameter space allowed by the observed properties of the newly 
discovered SNR GRO/RX J0852, assuming that the COMPTEL 44 Ti source and the ROSAT 
SNR are the same source. We conclude that the likely progenitor of this new SNR is a massive 
star and its precise type depends on the correct determination of the current SNR shell expansion 
speed v s . It is most probably that the true value of v s is more than a factor of 2 greater than that 
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derived from the shell electron temperature using X-ray data and the progenitor is a 15M SNII 
with a moderate SNE kinetic energy. The age of the source is ~ 700 yr and the distance is ~ 150 
pc. If the current v s were small, on the other hand, a much more energetic SNIb or SNIc event 
would be required to have exploded 800-900 yr ago in a high density medium about 250 pc away. 
This latter interpretation, however, has the difficulty of finding observationally the required high 
density gas at the required distance. The recently discovered cold neutral gas structure associated 
with the nearby Gum Nebula seems not to be able to provide the required high density. 

It is likely that the nearest SNE less than 1000 yr ago in the Vela region produced both the 
44 Ti and 26 Al emission features detected by COMPTEL. This offers an exciting new avenue to 
test the core-collapse SNE models since the explosive yields of 44 Ti and 26 Al are related. Because 
of its much shorter lifetime, 44 Ti is usually better suited for discovering young SNRs. However, 
since it just so happened that the age of GRO/RX J0852 is about 8-10 times the lifetime of 44 Ti, 
its 44 Ti flux has dropped by a factor of 3000-10,000 from its peak value and now almost equals 
that of 26 Al which has a decay lifetime that is ~ 10, 000 times longer. It is thus worth noting that 
for SNRs older than about 1200 yr, 26 Al is a better radionuclide to use for SNR searches. 

We thank Roland Diehl for communicating with us about the new analysis of the COMPTEL 
1.8 MeV feature in Vela and are very grateful to the anonymous referee for his/her insightful 
comments. 
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